Time-optimal synthesis of SU(2) transformations for a spin- 1/2 system 
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We consider a quantum control problem involving a spin- 1/2 particle in a magnetic field. The 
magnitude of the field is held constant, and the direction of the field, which is constrained to lie 
in the x — y plane, serves as a control parameter that can be varied to govern the evolution of the 
system. We analytically solve for the time dependence of the control parameter that will synthesize 
a given target SU (2) transformation in the least possible amount of time, and we show that the time- 
optimal solutions have a simple geometric interpretation in terms of the fiber bundle structure of 
SU(2). We also generalize our time-optimal solutions to a control problem that includes a constant 
bias field along the z axis, and to the case of inhomogeneous control, in which a single control 
parameter governs the evolution of an ensemble of spin-1/2 systems. 

PACS numbers: 42.50.Dv, 02.30.Yy, 02.20.Sv 



I. INTRODUCTION 

Many applications rely on the ability to coherently con- 
trol the state of a quantum system [lHa]- In particular, 
the current push to develop robust quantum information 
processors has led to the development of quantum control 
protocols for a diverse array of experimental platforms, 
including atomic, optical, and condensed matter systems 
• I n a typical control problem, the system in ques- 
tion is described by a Hamiltonian containing several con- 
trol parameters that we are free to vary, and we would 
like to determine the time dependence of these parame- 
ters such that the evolution of the system implements a 
desired unitary transformation. Such problems are gen- 
erally highly nontrivial: they do not usually admit an an- 
alytic solution, and must be solved via numerical searches 
[ll| - [l8| . Analytic solutions can, however, sometimes be 
obtained for control problems involving low-dimensional 
systems. In particular, for control problems involving a 
spin-1/2 particle, analytic solutions have been obtained 
that minimize either an energy-type cost functional [19l — 
[2~j| or the total evolution time [22[ HH . 

Here we consider a model quantum control problem 
involving a spin-1/2 particle in a magnetic field. The 
magnitude of the field is held constant, and its direction, 
which is constrained to lie in the x — y plane, serves as 
a control parameter that can be varied to govern the 
evolution of the system. The evolution can be described 
in terms of an SU{2) evolution operator U(t), such that 
if the state of the spin at time zero is 1^(0)) then the 
state at time t is \tp(t)) = U(t)\ifj(0)). Given an arbitrary 
target SU(2) transformation V, we analytically solve for 
the time dependence of the control parameter such that 
U(t) = V and t is as small as possible. By viewing SU (2) 
as a U(l) fiber bundle over the two-dimensional sphere 
S 2 , we are able to give a simple geometric interpretation 



to these time-optimal solutions. We also generalize our 
time-optimal solutions to a control problem that includes 
a constant bias field along the z axis. 

An important development in the field of quantum con- 
trol is the notion of inhomogeneous control, in which a 
single set of control parameters governs the evolution of 
an ensemble of systems subject to different Hamiltonians. 
The differences in the Hamiltonians may, for example, 
describe unwanted perturbations that give rise to deco- 
herence. By choosing the control parameters properly, 
one can compensate for these perturbations so that the 
result ing sy stem dynamics are insensitive to their pres- 
ence [2J, [25| . Alternatively, the differences in the Hamil- 
tonians may be intentional, so as to provide a means of 
addressing individual systems in the ensemble 26-28]. 

We investigate inhomogeneous control in our model 
control problem by generalizing the problem to the case 
of an ensemble of N spin-1/2 systems. The magnetic 
fields of the different systems vary in magnitude but are 
all aligned along a common direction in the x — y plane, 
and we take this common direction to be the control pa- 
rameter that governs the evolution of the entire ensemble. 
We obtain a semi-analytic solution to this inhomogeneous 
control problem for the case N = 2, and we verify that 
our solution is time-optimal by comparing it with the 
results of a numerical search. 



II. CONTROL PROBLEM 

The system that we consider consists of a spin-1/2 par- 
ticle in a magnetic field B. We assume that the magni- 
tude B = \B\ of the magnetic field is constant, and its 
direction n = B/\B\ serves as a control parameter that 
can be varied to govern the evolution of the system. The 
Hamiltonian for the system is 



H = —fiBer ■ n, 



(1) 
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where n is the magnetic moment of the particle and at 
are the Pauli spin matrices. For simplicity, we will choose 
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units such that /j,B — 1. The system evolves in time 
according to the unitary transformation 



U(t)=Texp(-i [ H{t')dt'), 
Jo 



(2) 



where T is a time-ordering operator that places operators 
at early times to the right of operators at later times. We 
note that U satisfies the Schrodinger equation 



iU = HU. 



(3) 



From Eq. ©, and the fact that H is traceless, it follows 
that det U = 1, so U is an SU(2) transformation. 

We now consider a control problem in which we are 
given a target SU(2) transformation V and are asked to 
determine the time dependence of the control parameter 
n and the total evolution time t such that U(t) = V 
and t is as small as possible. If h is allowed to point in 
any direction, then the solution to the control problem is 
trivial: we write V in the form V = e lr ' CT , where \r\ < tt, 
and we take 



n = r, 



t 



(4) 



For example, for a target transformation V = e" 7 "'*' 2 
describing a spatial rotation with axis z and angle r], we 
find that n = z and t = rj/2. 

Let us suppose, however, that the control parameter 
n is constrained to lie in the x — y plane. The control 
problem is still solvable, but the solution is no longer 
trivial. We can verify that the control problem is solvable 
by presenting a solution that is not time-optimal. Let 
us write the target transformation V in terms of Euler 
angles 0, and <j>: 



.j/2 i<f>cr x /2 



(5) 



From Eq. ©, it follows that V can be synthesized by by 
taking 



n(r) 



for < r < \<f>/2\, 

for |0/2| < r < |0/2| + \e/2\, (6) 
for 10/21 + \6/2\ <T<t, 

t=\1>/2\ + \6/2\ + \<t>/2\. (7) 



y 

X 



For example, consider again a target transformation V = 
gMjtfz /2 describing a spatial rotation with axis z and angle 
77. We find that V = e -*W V^/V™*/ 4 , so = -ip = 
7r/2, 9 = rj, and t = ir/2 + r]/2. For comparison, recall 
that t = 77/2 for the unconstrained control problem in 
which n is allowed to point in any direction. 



III. TIME-OPTIMAL SOLUTION 

We now present a time-optimal solution to the con- 
strained control problem. We begin by describing two 
methods for assigning coordinates to an SU (2) transfor- 
mation U. For the first method, we assign real-valued 



coordinates r = (w,x,y, z) to U by expanding U in the 
Pauli spin matrices: 



U = w + ixa x + iyo y + iza z . 



(8) 



We call these coordinates embedding coordinates, be- 
cause they describe an embedding of SU(2) into R 4 . 
For the second method, we assign complex-valued coor- 
dinates (#1,122) to U by expressing U in the form 



U 



zi z 2 



(9) 



We call these coordinates complex coordinates. From 
Eqs. © and ©, it follows that the two sets of coordinates 
are related by (z±, Z2) = (to + iz.y + ix). 

The Lie group SU(2) is three-dimensional, but both 
sets of coordinates label SU(2) transformations using 
four real parameters. So for both sets of coordinates 
there are more coordinate degrees of freedom than phys- 
ical degrees of freedom, and only some of the points in 
the coordinate space actually correspond to SU(2) trans- 
formations. From Eqs. © and ©, it follows that such 
points satisfy the constraint 



Z2 



1. 



(10) 



The locus of points r that satisfy Eq. (flQ|) is a three- 
dimensional sphere S 3 embedded in 1R , and the mapping 
U h-> r is a diffeomorphism from SU{2) to S 3 . 

It is useful to express the Schrodinger equation © 
in terms of both sets of coordinates. We first consider 
the embedding coordinates. We substitute the definition 
of the embedding coordinates given in Eq. © into the 
Schrodinger equation © to obtain an equation of motion 
for r: 



r = n x L x (r) + n y L y (r) +n z L z (r) 7 (11) 



where 



L x (r) = wx — xw + zy — yz, (12) 
L y (r) = wy — yw + xz — zx, (13) 
L z (r) = wz — zw + yx — xy, (14) 

are orthonormal vectors that span the tangent space of 
S 3 at the point r. As r evolves in time, it traces out a 
path in S* 3 whose tangent vector is r. From Eq. ([TT]) it 
follows that the length of the tangent vector is \r\ = 1, 
so time corresponds to arc length along the path. The 
time evolution of r is governed by the control parameter 
n, which dictates the projection of the tangent vector r 
along the basis vectors Lf-(r): 



L x (r) ■ r = 
L y {r) ■ r = 
L z (r) ■ r = 



(15) 
(16) 
(17) 
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For the constrained control problem h = cos 4>x + sin <fi y 
for some angle </>, so 

L x (r) ■ f = cos 4>, (18) 
Ly(r)-r = aw<l>, (19) 
L»(r)-r = 0. (20) 

It is also useful to express the Schrodinger equation 
§3§ in terms of the complex coordinates. We substitute 
the definition of the complex coordinates given in Eq. ([9]) 
into Eq. ([3]) to obtain equations of motion for Z\ and Z2'- 

zi = -ie'^ z%, (21) 
Z2 = ie-^zt. (22) 

If we differentiate Eqs. (f2"Tj) and (|2"21 with respect to t and 
then substitute for z\ and 22 using the original equations, 
we obtain the decoupled equations 

'i\ + itfizi + z\ = 0, (23) 
z 2 + + z 2 = 0. (24) 

Using Eqs. ([2"T]) and ([2"2"]h it is straightforward to derive 
the identities 

zjzj + z 2 z* 2 = 0, (25) 
ziz x - iiz 2 = ic - **. (26) 

We can understand the meaning of these identities by 
transforming from complex coordinates to embedding co- 
ordinates: 

z\z\ + i 2 z 2 — r ■ r + iL z (r) ■ r, (27) 
z 2 z 1 - z x z 2 = L y (r) ■ r + iL x (r) ■ r. (28) 

So Eqs. and @ follow from Eqs. @DJ) and ([T5 ]) -((5D ]> . 

Let us now return to the embedding coordinates and 
consider the problem of finding a minimum-length path 
in S 3 that satisfies the constraint L z (r) ■ r = 0. Such a 
path can be obtained by minimizing the action 

S = J (\r'\ + 7 (|rf - 1) + XL z (r) ■ r') du. (29) 

Here u is an arbitrary parameterization of the path, 
r' = dr/du, and 7 and A are Lagrange multipliers. The 
first term of the integrand gives the length of the path, 
the second term imposes the constraint |r| 2 = 1, which 
restricts the path to S 3 , and the third term imposes the 
constraint L z (r)-r — 0, which expresses the fact that the 
control parameter n must lie in the x — y plane. Note that 
r' = dr/du = (dt/du)r and \r\ = 1, so the parameter u 
is related to the time t by 

dt/du=\r'\. (30) 

We write down the Euler-Lagrange equations correspond- 
ing to the action given in Eq. (|2TJ)) . use Eq. (|3H)) to replace 



u with f, and transform from embedding coordinates to 
complex coordinates to obtain 

z\ + 2i\zi + (iX - 27)21 = 0, (31) 
z 2 + 2i\z 2 + (iX - 2-f)z 2 = 0. (32) 

A time-optimal solution to the constrained control 
problem must satisfy the Schrodinger equations (|2"3"|) and 
as well as the Euler-Lagrange equations (|3"Tj) and 
([32]). We subtract Eq. ([23]) from ^ and Eq. ([24]) from 
(1521 to obtain 

i(2X - 4>)zi + (iX - 27 - l)zi = 0, (33) 
i(2X - 4>)z 2 + (iX - 2 7 - l)z 2 = 0. (34) 

Using the identities given in Eqs. ([23]) and ([2"o) , we can 
eliminate the coordinates Z\ and z 2 from Eq. ([33]) and 
(15^1) and obtain equations that involve only the parame- 
ters 7, A, and 0: 

= 2A, iA = 2 7 + l. (35) 

The solution to these equations is 

7 = -1/2, A = w/2, = <f> + ujt, (36) 

where 4>q and w are integration constants. So an SU (2) 
transformation can be synthesized in a time-optimal fash- 
ion by varying the control parameter <fi as described by 
Eq. ([36]). 

We would now like to calculate the evolution operator 
U that results when the control parameter <p is varied 
in the time-optimal fashion described by Eq. ([33]) . We 
first note that t/(0) is the identity transformation, which 
has complex coordinates (2:1,2:2) = (1,0). We substitute 
Eq. ([33]) for tfi into the Schrodinger equations (|2"TT) and 
([2"2"]l and solve them subject to these initial conditions to 
obtain 

21 = (2a)- 1 ((3 + e^- t + ^e" 4 ^*), (37) 

2 2 = (2ay 1 e' l4 ' (e^- t - e" l/3+ *), (38) 

where 

a = (1 +w 2 /4) 1/2 , P±=ol±uj/2. (39) 

It is useful to view the parameters (</>o, w, i) as defining a 
third set of coordinates for U . We call these coordinates 
time-optimal coordinates. Eqs. (|37p and ([55]) can then 
be viewed as describing a coordinate transformation from 
time-optimal coordinates to complex coordinates. 

Suppose we are given a target SU{2) transformation 
V. We can synthesize V in a time-optimal fashion by 
determining its complex coordinates (21, 22) and then in- 
verting Eqs. ([37]) and (]38l) to obtain its time-optimal co- 
ordinates (</)o,w,t). The parameters (fio and ui tell us the 
time dependence of the control parameter </>, and the pa- 
rameter t tells us the total evolution time. 
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FIG. 1: (Color online) Time t needed to synthesize the trans- 
formation V = e l * ,<Tz ' /2 versus r\. (a) Time-optimal solution for 
the constrained control problem, (b) Euler solution for the 
constrained control problem, (c) Time-optimal solution for 
the unconstrained control problem. 



Let us now consider some specific examples. First we 
consider a target transformation V = e"' ee ' CT / 2 that de- 
scribes a spatial rotation with axis eg = cos 6 x + sin 8 y 
and angle 77. The complex coordinates of V are (z\, z 2 ) — 
(cos r//2,ie- ie sin ry/2). We invert Eqs. ([37]) and ^ to 
obtain the time-optimal coordinates: 



u> = 0, 



t = v /2. (40) 



This solution is identical to the time-optimal solution for 
the unconstrained control problem described in Eq. (|4|). 
This is to be expected, since the time-optimal solution for 
the unconstrained control problem satisfies the constraint 
that n must lie in the x — y plane. 



Next we consider a target transformation V — e 



irj(T z /2 

that describes a spatial rotation with axis z and angle 
rj. The complex coordinates of V are (z i7 Z2) — (e J?? / 2 , 0). 
We invert Eqs. (J37J and ([55]) to obtain the time-optimal 
coordinates: 



= 2 V {l-v 2 )- 1 '\ t = n(l-v 2 ) 1 / 2 , 



(41) 



where v = 1 — 77/271". The parameter tfrg is undetermined 
by the inversion, and any value can be used to perform 
a time-optimal synthesis of V. Mathematically, 4>q is un- 
determined because V is located at a coordinate singu- 
larity of the time-optimal coordinate system; physically, 
it is because V is invariant under similarity transforma- 
tions involving arbitrary rotations about the z axis. In 
Fig. [T] we compare the time-optimal solution described in 
Eq. (|4"Tj) with the Euler solution described in Eqs. 
(UJ and the time-optimal solution for the unconstrained 
control problem described in Eq. (U]). 

Let us now consider the trajectory of the spin on the 
Bloch sphere as it evolves along a time-optimal path. If 
the state of the spin at time zero is \ip(0)), then the state 
at time t is \ip(t)) — U(t)\ip(0)) . We can represent the 



state of the spin at time rj as a point s(t) = (ip(t)\a\i[i(t)) 
on the Bloch sphere. In Fig. [2] we plot the trajectory of 
the spin on the Bloch sphere for the time-optimal syn- 
thesis of a 7r/2-rotation about the 2 axis 
where the spin is initially aligned along the i axis for 
Fig-d^a) and the — y axis for Fig.[5Jb). For both curves 
we take <po = 0. 



IV. PROPERTIES OF THE TIME-OPTIMAL 
SOLUTIONS 

We can visualize the time-optimal solutions by rep- 
resenting SU(2) transformations as points on the two- 
dimensional sphere S 2 . Given an SU(2) transformation 
U, we define p(U) to be the point on S 2 corresponding 
to the state W \ t) : 



p(U) = <t \U*rf\ t). 



(42) 



We note that p(U) = p{e l8,T -U) for any value of 6. This 
property of p allows us to view SU(2) as a fiber bundle, 
where S 2 is the base manifold, U(l) is the fiber, and 
p : SU(2) — > S 2 is the projection function. 

We will now show that the time-optimal solutions 
project to circles on S 2 . Let us identify the plane that 
bisects S 2 at the equator with the complex plane. We 
can map points p on S 2 to complex numbers C(p) 011 the 
complex plane by stereographically projecting from the 
south pole: 



C(p) = 



Px + iPy 
1+Pz 



(43) 



Let (zi,Z2) denote the complex coordinates of an arbi- 
trary SU(2) transformation U. From Eqs. ©, (|4"2"]l. and 
(1331. it follows that 



C(p(U)) = Zl /z 2 . 



(44) 



For a time-optimal solution, z\ and z 2 are given by 
Eqs. (|3"Tj) and (|3"5)l . We substitute these expressions into 
Eq. g4]) to obtain Q(t) = f(e 2lat ), where 



(3+z + f3- 



(45) 
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FIG. 3: (Color online) Paths on the two-dimensional sphere 
S 2 for the time-optimal synthesis of the transformation V = 
e rr > a z' 2 ) where rj = n/2, n, 3n/2, 2n. Longer paths correspond 
to larger values of r\. 

The function f(z) is a Mobius transformation. Since e 2lat 
describes a circle in the complex plane, and both stereo- 
graphic projection and Mobius transformations preserve 
circles, it follows that the time-optimal solutions project 
to circular paths on S 2 . In Fig. [3] we plot example 
paths for the time-optimal synthesis of the transforma- 
tion V = e Lr, ' Tz / 2 , which describes a spatial rotation with 
axis z and angle ij. The paths begin and end at the north 
pole. For the paths shown we take fa = 0; alternative 
paths that also synthesize V can be obtained by taking 
different values of fa, and for such paths Fig. [3] is rotated 
about the z axis through an angle fa. Under the fiber 
bundle interpretation, the time-optimal solutions can be 
obtained by lifting the circular paths from S 2 to SU(2), 
where the lifts are performed relative to the connection 
induced by the constraint L x (r) ■ r = 0. Another way 
to visualize the time-optimal solutions is to stereographi- 
cally project from the north pole, in which case the time- 
optimal solutions map to straight lines on the complex 
plane. 

We have shown that time-optimal solutions project to 
circular paths on S 2 . We will now show that the length 
of the path on S 2 is equal to twice the amount of time 
needed to synthesize the corresponding transformation. 
We first assign coordinates (ip, 9, fa) to an arbitrary SU (2) 
transformation U by performing an Euler-angle decom- 
position: 

U — e iV<r 2 /2 e i0<T H /2 e i0<T 2 /2^ /^g\ 

We call these coordinates Euler coordinates. Note that 

p(U) = sin 9 cos fax + sin 9 sin fay + cos 9 z, (47) 

so the coordinates (9, fa are the spherical-polar coordi- 
nates of the point p(U) on S 2 . From Eqs. flSJ), © and 
(|46[). it follows that the Euler coordinates are related to 
the complex coordinates (zi,z 2 ) and the embedding co- 
ordinates r = (w, x, y, z) by 

z 1 =w + iz = e l ^ +<t ' )/2 cos9/2, (48) 

z 2 = y + ix = e l( ^ )/2 sin(9/2. (49) 

Let us consider a small segment [t, t+dt] of a time-optimal 
path on S 3 . From Eq. (|48p. it follows that the arc length 



dt of the segment is given by 

dt = (dr ■ dr) 1/2 

= (l/2)(d9 2 + d(f> 2 +dijj 2 + 2cos9dfaifa 1/2 . (50) 

Recall that time-optimal paths satisfy the constraint 
L z (r) r = 0. From Eqs. an d (0H|j it follows that in 
Euler coordinates this constraint takes the form 

dip + cos9d<p = 0. (51) 

We substitute Eq. (|5TJ into Eq. to obtain 

dt 2 = (1 / A)(d9 2 + sin 2 9 dfa) = (1/4) ds 2 , (52) 

where ds 2 is the standard metric on S 2 , which is induced 
by the Euclidean metric on R 3 via the embedding of 5* 2 
into R 3 . From Eq. (|52l) . it follows that the time needed 
to synthesize an SU(2) transformation is equal to half 
the length of the corresponding path in S 2 . 

V. BIAS FIELD 

Let us now generalize the control problem described in 
Sec. [IT] by adding a constant bias magnetic field along the 
z axis. The Hamiltonian for the system is now given by 

H = cr ■ n + ba z , (53) 

where b characterizes the strength of the bias field. As 
before, we assume that n is constrained to lie in the x — y 
plane and thus has the form h = cos fax + sin. fay. We 
assume that we are given a target SU(2) transformation 
V and bias field value 6, and we would like to determine 
the time dependence of <f> and total evolution time t so 
as to synthesize V in a time-optimal fashion. 

It is convenient to work in the interaction picture. We 
express the Hamiltonian as H — Hq + Hi, where Ho = 
ba z is the bare Hamiltonian and Hi = cr ■ n is the 
interaction Hamiltonian, and we define Ui = e lHot U to be 
the interaction-picture evolution operator. The operator 
Ui satisfies the Schrodinger equation 

ill, - HrUi, (54) 

where 

Hi = e iHot Hie- iHot = -a ■ hi, (55) 

hi = x cos fa- + y sin fa , (56) 

fa = 4>+ 2bt. (57) 

From the results of Sec. IIII1 it follows that the time- 
optimal solution for fa- is given by fa- — fa + ujt, where 
fa and ui are constants, and the complex coordinates 
(z 1 (Ui),z 2 (Ui)) of fa are given by Eqs. ([37]) and ([38]>. 
Since U — e~ lHot Ui, it follows that the complex coordi- 
nates (zi(U), z 2 (U)) of U are given by 

Zl (U) = (2a)- 1 e- M (/3 + e lfl - t + ^e^ 4 ), (58) 
z 2 {U) = {2a)- 1 e-' l ^ o+bt \e l0 - t - e~ l0+t ), (59) 
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where a and f}± are given by Eq. (|39|) . Given the complex 
coordinates of the target transformation V, we can invert 
Eqs. ([58]) and ([59]) to determine the parameters needed 
to synthesize V in a time-optimal fashion. 



VI. INHOMOGENEOUS CONTROL 

We will now generalize the control problem described 
in Sec. [II] to the case of inhomogeneous control. We con- 
sider an ensemble of N spin-1/2 particles, where particle 
i is in a magnetic field B L = Bih with magnitude Bi and 
direction n. The Hamiltonian for particle i is 



Hi = 



(60) 



where Xi = I^Bi. As before, we assume that n is con- 
strained to lie in the x — y plane and thus has the form 
n = cos (f>x + sin <j> y. We note that the single control pa- 
rameter cf> governs the evolution of all N particles. If we 
evolve the ensemble for a time t while varying the con- 
trol parameter <fi, we obtain SU (2) evolution operators 
{Ui(t),--- ,C/Ar(t)}, where Ui(t) is the evolution opera- 
tor for particle i. We assume that we are given a list of 
target SU(2) transformations {Vi, • • • , Vjv} and a list of 
field values {xi, ■ • ■ ,Xn}- We would like to determine 
the time dependence of (f> and total evolution time t such 
that Ui(t) = Vi for i = 1, • • ■ ,N, and t is as small as 
possible. 

We begin by adapting the formalism developed in 
Sec. IHII to the case of the Hamiltonian Hi given in 
Eq. ([6"0"]) . We denote the embedding coordinates of U, 
by Ti and the complex coordinates of Ui by (zu, Z2i)- 
The Schrodinger equation in embedding coordinates is 



Xi{n x L x (r.i) + n y L y (rA + n z L z (ri)), 



(61) 

From Eq. (|6ip and the orthonormality of the vector fields 
Lk, it follows that the magnitude of the tangent vector 
Ti is | ri| = Xi, so the arc length s of the path traced 
out by Ti in S 3 is related to the time t by s — \it- The 
Schrodinger equation in complex coordinates is 

(62) 
(63) 

From Eqs. (|6"3"]) and (f6"3"|) we obtain the decoupled equa- 
tions of motion 



Zli 



'in + i<t>hi + XiZu = 0, 
z 2i + i(j>z 2 i + x\z-2i = 0, 
and the identities 



Z\iZ* u + Z 2l Z2i = 0, 

Z2iZu - z u z 2i = iXid 



-i(j> 



(64) 
(65) 

(66) 
(67) 



We can obtain a time-optimal solution to the control 
problem by minimizing the action 



S 



(68) 



where 

Ai=Xil (HI + 7i(W 2 - 1) + KL z (n) ■ <) du, (69) 



J (L x (ri) ■ r[ - L x {r.j) ■ r\) du, 
C tj = Cij J {L y (n) ■ r'i - L y ( rj ) ■ r'j) du, 



(70) 
(71) 



and 7,, Xi, bij, and c,j are Lagrange multipliers. The 
terms Ai are straightforward generalizations of the action 
(|29"1) for the original control problem; the prefactor Xi 
accounts for the fact that the arc length s of a path in 
S 3 is related to the time t by s — xrf- The terms By 
and Cij impose the constraints L x (r.i) ■ Vi = L x (rj) ■ Vj 
and L y (ri) ■ r'i — L y {r.j) ■ rf, from Eqs. (fT5j) and (TT5]), we 
see that these constraints account for the fact that the 
same control parameter (f> governs the evolution of all N 
evolution operators {U\, . . . ,Un} ■ 

We now follow the same procedure described in 
Sec. IIIII we write down the Euler-Lagrange equations 
for S, subtract the decoupled Schrodinger equation (|6"2")l 
and (|63|) . and use the identities (|66l) and (|67|) to obtain 
equations that involve only the Lagrange multipliers and 
the control parameter </>. We find that 



x: 



'■(2X l -<j } )^e^J2(w lJ ~w ji ), (72) 



i\, - 2 7l - X 2 i = -2ie** J2(w tJ - Wji ), (73) 



where w%j = bij + ifyj . 

For the case N = 2 we can solve Eqs. (J72J) and ([731) 
to obtain an equation of motion for (j>. From Eqs. ([721) it 
follows that 



Ai = (l/2)W>- 
A 2 - (l/2)(0 



oc/xl), 
a/xl), 



(74) 
(75) 



where 



J4> 



a = we"^ (76) 

and w = W12 — tf2i- From Eqs. (|73[) it follows that 

Ai + A 2 = 0, (77) 
u; = -(l/4)(Ai-Aa + 2i/3)e- i * (78) 

where /3 = 2^ x +x\ = ~{2^ 2 + xl)- We integrate Eq. (J77J 
to obtain 



Ai + A 2 = A, 



(79) 



where A is an integration constant. We solve Eqs. (|74l) . 
([75]). and ([751) f° r Ai, X 2 , and a in terms of 4> an d A: 



x 1 = ( x /2)<Mx72x?)A 

A 2 = -(x/2)0+(x72xl)A 



(80) 
(81) 
(82) 
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where 



2 i 2 
Xl+X2 

— 9 05 

Xl -%2 



x ' ee (83) 
X1-X2 



We substitute Eqs. ([50)1 and (I5T1) for Ai and A 2 into 
Eq. (J7HJ) to obtain 



w = -(l/4)(x<?i + 2i/3)e-^. 



(84) 



We differentiate Eq. with respect to time and sub- 
stitute the resulting expression for w into Eq. (|76[) to 
obtain 



a = -(1/4)( X + 2i(3 - + 20)). 



(85) 



Taking the real and imaginary parts of Eq. (|85j) , we find 
that 



(86) 
(87) 



a = -(l/4)(x0+2/30), 
= -(l/4)(2/3-x#). 
We integrate Eq. (|87|) to obtain 

= (x/4)0 2 + B, 



where B is an integration constant. Substituting 
Eqs. for a and §53) for /3 into Eq. (JHSJ) , we find 
that 

"4> + 3 /2 + (2B/x)0 + (4x'/x)(0 - A) = 0. (89) 

So the control parameter satisfies the equation of mo- 
tion 



(j) +<f> 3 /2 + b<]) + a = Q, 



(90) 



where a = -(4x'/xM and 6 = (2/x)S-4%'/X- We note 
that since the integration constants A and B can take 
any values, the parameters a and b can also take any 
values, and are thus not constrained by the values of xi 
and X2- 

Given initial conditions (0o, </>o, 4>o) and parameters 
(a, b), we can integrate Eq. (|90p to obtain a time-optimal 
solution for </>. Given this time-optimal solution, we can 
integrate the Schrodinger equations (|64| and (|65l) sub- 
ject to the initial conditions (zu,Z2i) = (1,0) to ob- 
tain the complex coordinates of a pair of evolution op- 
erators {£/i,C/ 2 }. It is useful to view the parameters 
(</>o, <fia, </>o, a, b, t) as a generalization of the time-optimal 
coordinates described in Sec. Mil The two integrations 
then define a coordinate transformation from the time- 
optimal coordinates to the complex coordinates of the 
pair of evolution operators {U±, 1/2}- Given target SU (2) 
transformations {V±, V2} and field values {xi, X2}, we can 
write down the complex coordinates of {Vi, V2} and then 
invert this coordinate transformation to determine the 
time dependence of the control parameter </> and the to- 
tal evolution time t needed to synthesize V± and V2 in a 
time-optimal fashion. We have thus formally solved the 
inhomogeneous control problem for the case N = 2. 




FIG. 4: (Color online) Control parameter <j> versus time 
t. The solid curve is obtained by numerically integrating 
Eq. (|90[) ; the points are obtained from a numerical gradient- 
ascent search with t — 3 and R = 50. 



We note that the parameters (</>o, 4>o, 4>o, a, b) deter- 
mine a time-optimal evolution for the control param- 
eter </>, and this evolution, together with the parame- 
ters (t, xi, X2)j determine a pair of evolution operators 
{1/1,1/2}. It is interesting that the time-optimality of 4> 
does not depend on the field values xi and X2- That is, 
if we hold the time dependence of <f> and the total evolu- 
tion time t fixed, and vary \i and X2, we will synthesize 
different evolution operators U\ and U2, but it will al- 
ways be the case that the synthesis of these operators is 
time-optimal. 

Let us now consider a specific example. We will take 
the field values to be xi = 1/2 and X2 = 3/2, and con- 
sider the pair of transformations {Vi, V2} whose time op- 
timal coordinates are </>o = 0, 4>q = —2, tpo = 0, a = 2, 
b = 3, t = 3. We numerically integrate the equation 
of motion ([90]) to determine the time evolution of the 
control parameter <fi that synthesizes Vi = U±{t) and 
V2 = L^2 (t) in a time-optimal fashion, and we numerically 
integrate the Schrodinger equations (|rJ2"j) and (|6"3")l to de- 
termine the complex coordinates of the pair {V±, V2}. In 
Fig. H] we plot the resulting time-optimal evolution of <fi. 

We verify that the synthesis of V\ and V2 is time op- 
timal as follows. Given arbitrary SU{2) transformations 
A\ and A2, we define the fidelity with which A\ and A2 
approximate V\ and V2 to be 



(l/A)(Tr[V^A 1 ] 



Tt[V]A 2 ] 



(91) 



The fidelity ranges from —1 to 1, where F = 1 if A\ = Vi 
and A2 = V2, and F decreases as the deviation of A\ 
and A2 from V\ and V2 increases. We fix the total evo- 
lution time i, and we discretize the time evolution of the 
control parameter by dividing t into R timesteps of dura- 
tion St = t/R. We define <j) r = cf){r8t) to be the value of 
the control field at timestep r. We then take A\ = U\(t) 
and A2 = U2 (t) , and perform a numerical gradient-ascent 
search to maximize F with respect to the discretized con- 
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FIG. 5: (Color online) Maximum fidelity F max versus time 
t, as determined by a gradient-ascent search with R = 50. 



tion of <fi for t = 3, as determined by the gradient-ascent 
search. We find good agreement with the time-optimal 
evolution of <j> obtained by integrating the equation of 
motion (JSDJ- 



VII. SUMMARY 

We have considered a quantum control problem in- 
volving a spin-1/2 particle in a magnetic field. We have 
analytically solved for the time dependence of the con- 
trol parameter needed to synthesize an arbitrary SU{2) 
transformation in a time-optimal fashion, and we have 
generalized our solution to the case of an inhomogeneous 
control problem involving an ensemble of spin-1/2 sys- 
tems. 



trol parameter values {<j>Q, ■ ■ ■ , 4>b~i}- In Fig. [S] we plot 
the numerically-determined maximum fidelity F max as a 
function of t for R = 50. Since F max first reaches 1 at 
t = 3, we see that the evolution described above is indeed 
time-optimal. In Fig. [H we plot the time-optimal evolu- 
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